'Electron' and 'photon' emerging from supersymmetric neutral particles: A possible 
realization in ultracold Bose-Fermi atom mixture 
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We show that the 'electron' and 'photon' can emerge from a supersymmetric Hubbard model 
which is a non-relativistic theory of the neutral particles. The Higgs boson and 'photon' may not 
appear in the same phase of the phase diagram. In a Mott insulator phase of the boson, the 
'electron' and 'photon' are stablized by an induced Coulomb interaction between 'electrons'. This 
emergent mechanism may be 'realized' in an ultracold Bose-Fermi atom mixture except the long 
range Coulomb interaction is repalced by a nearest neighbor one. We suggest to create 'external 
electric field' so that the 'electron' excitation can be observed by measuring the linear density-density 
response of the 'electron' gas to the 'external field' in the time flying experiment of the mixture. 
The Fermi surface of the 'electron' gas may also be expected to be observed in the time flying. 

PACS numbers: 03.75.Lm,12.90.+b,11.15.-q,71.30.+h 



What roles did the mysterious Higgs boson play ex- 
cept providing a mass to the intermediate bosons in the 
weak-electrodynamics? Where did the electron and pho- 
ton come from [lj? In this Letter, we suggest that it was 
a neutral world in high temperature and the electrons 
and photons emerged when the temperature was lower 
than certain critical temperature. The Higgs boson fluc- 
tuation plays an important role in the appearance of the 
finite gauge coupling constant. We show that this emer- 
gence can be 'realized' by the elementary excitations in 
the ultracold Bose-Fermi atom mixture in optical lattices 
except the Coulomb interaction between 'electrons' is re- 
placed by a nearest neighbor one. 

We start from a three-dimensional spinless supersym- 
metric(SUSY) Hubbard model, which was called 'ultra- 
cold superstring' in one-dimensional lattice Q. The mi- 
croscopic description of the model parameters for an ul- 
tracold atom mixture has been proposed [j| . The consti- 
tution particles in this model are neutral spinless bosons 
ifli) and fcrmions (/j) with i the lattice site index. In 
high temperatures, there are only excitations of these 
constitution particles. We call this a confinement phase. 
After the SUSY is broken, the system undergoes a phase 
transition, in certain critical temperature, to a uniform 
mean field (UMF) state in which a finite coupling gauge 
field and a 'charged' fermion emerge if the fermion occu- 
pation is nearly half filling. We call these emergent ob- 
jects 'photon' and 'electron', respectively. At the exact 
half filling, this UMF state turns to a long range ordered 
state, the checkerboard crystal 0. However, this UMF 
state can only be stable if the lattice filling of the con- 
stitution bosons is integer and has a Mott insulator (MI) 
ground state. For a special space structure, the 'electron' 
interaction is pure Coulomb. 

Experimentally, the Bose-Fermi atom mixture in opti- 
cal lattices has been realized for 87 Rb- 40 K and 23 Na- 
6 Li 6]. We suggest an experiment to create an 'external 
field' by changing the depth of the fermion's optical po- 
tential. The response function of the mixture to the ex- 



ternal field may be measured by the density distribution 
in a time flying of the cold atom cloud. The behavior 
of the response functions may be used to identify the 
fcrmionic elementary excitations in a set of given param- 
eters, either the constitution one or 'electronic'. On the 
other hand, we expect the 'electron' Fermi surface can 
be observed by the time flying, which has been used to 
observe the Fermi surface of the pure cold Fermi atoms 


The Bose-Fermi mixture Hubbard Hamiltonian we are 
interested in reads 



H = - //;.,,«;«., + t F ,ijf-fj) - ^(a«b< + fJ.Fi 
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where the lattice spacing is set to unit, nf = a\ai and 
n{ — fjfi- tB,ij and tp.ij are the hopping amplitudes of 
the boson and fermion between a pair of sites. [Ib and /if 
are chemical potentials. And Ubb and Ubf are the on- 
site interactions between bosons, and between boson and 
fermion, both of which may be repulsive or attractive. 
The microscopic calculations of these model parameters 
in terms of the cold atom mixture may be found, e.g, in 



If t B , 



t F A 



kj, fJB = ftp = fJ and Ubb 



Ref. 

Ubf = U, the model is SUSY under transformations 
Sai = iQfi, Sfi = iOa.i where 9 is a Grassmann number. 

To deduce the low energy theory in strong couplings , 
we will use the slave particle technique, which has been 
applied to the cold boson system |3, llCM. Ill| . In the slave 
particle language, the Hamiltonian reads H = H2 + H4 
where Hi and H4 are the two-operator and four-operator 
terms, respectively. Namely, 



(2) 
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and 

ij a, 



Aj w A (which is a saddle point value of A;). This is cor- 
responding to a mean field approximation. The effective 
mean field action is given by 



- E^iE^+V^+i (3) 

ij a, (5 
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where 



and 



l hi = h a .jha,i- We explain 



briefly the deduction of this slave particle Hamiltonian. 
The state configurations at an arbitrary given site con- 
sists of {\a, s) | a = 0, 1, 2, ... ; s = 0, f } where a and s 
are the boson and fermion occupations, respectively. The 
Bose and Fermi creation operators can be expressed as 
« f = E o V«+T[|a + l,0)<a,0| + |a + l,l><o,l|],/ t = 

Ea I". l>(a.0|. The mapping to the slave particle reads The mean fieJd equations are = 



where (D^r. 1 = (<9 r - n F )S tj - t F Y,fXa,ji 6 j,i+r, Mf = 
(iBa+VF + !^?r-ct(a-l) + UBFa~i\ and (-D^)^ 1 = (9 T - 

a^B + % g -a(a- 1) -iA))% -tp E-? Xoj^j.i+f- (t is the 
unit vector of the lattice.) To simplify, we have assumed 
only the nearest neighbor hopping tp of the fermions is 
not zero. 



,0) - hl,\aA) 



We call c a the composite 



"a,i c a,j) 



fermion (CF) |J] and h a the slave boson. The normal- 
ized condition E Q (l a ) 0| + |a, l)(a, 1|) = 1 implies a 
constraint E Q ( n c + n /J) = 1 at each site. The slave parti- 
cles arise a U(l) gauge symmetry c a , 



e l ' Pi h , i i. The SUSY transformation is given by Sh a ,i — 
i8y/ac a -i^, 5c a ,i — iQ\Ja. + lh a+ i^. The global U(l) 
symmetry c a ^ — ► e latp c ay i, h a ^ — > e %aif, h a ^ reflects the 
particle number conservation. Since the slave particle 
technique essentially works in the strong coupling region, 
we focus on Ubb = Ubb/\sq\ ^> 1 where Ek is the dis- 
persion of the constitution boson. (Hereafter, a quantity 
A = A/\e Q \.) 

There are two types of four slave particle terms in 
H4, the ts terms and tp terms. We first neglect the 
ts terms in the mean field level of the CF. To decou- 
ple the tp terms , we introduce Hubbard-Stratonovich 
fields tfufi ij an d Xa/3,ij- The partition function is given 
by Z = J DxDfjDcDcDhDhDX e~ S 'ff where the effec- 
tive action reads 



£e//[X,17,C,/l,A] 
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+Va0,ij(Xp a ji ~ hpjh a ,i)] + t B terms \, 



(4) 



where \ is the Lagrange multiplier for the constraint 
ScK. + n h,i) = L Rewriting = V^ tij e iAij and 

Xa0,ij = Xa0,ij eiAii and A, = A + Ao,i, Aoi and A i3 are 
U(l) gauge field corresponding to the gauge symmetry. 
Fixing a gauge, » V^Sap, Xap,ij ~ Xa,ij S <*P and 



T En D c i j(p n ), Xa, 



where _ £>£ tf (p„) = J dre*» T D^(r) and D^K) = 
/ dre lw " T D^ ij (T) with p n and w„ the Fermi and Bose 
frequencies. Solving the mean field equations near the 
critical point with \ a w 0, one can obtain the critical 
temperature T C Q 



(6) 



with Tp B 



TpXaji and T« X ^ 



n Q _ ^^ e Pl-i\-utiB+a(u-l)U B B /2] _ j^j n Q _ 
X^ e /3[-«A-(aMB+PF)+a(a-l)L r Bi3/2+Q(7BF] _|_ J]. Below 

this temperature, the minimized free energy including 



2(S 2 ° 



where t q — 



T?-T 



variables x is F cx — Eq( t ") 
and S 1 ^, are the fractions of non-zero terms corre- 
sponding to the second order and four order of \x a \ H|. 



The optimal \Xa\ 



"S%/S%. 



Before going to a concrete mean field state, we first re- 
quire the mixture is stable against the Bose-Fermi phase 
separation. For example, it was known that the mixture 
is stable in the Ml phase if Antp sm(-Kn^)UBB > U BF 0. 
For the parameters in the SUSY model, this inequality 
is not satisfied. Thus, we assume the SUSY has been 
broken. In cold atom contents, the SUSY is broken by 
hands through changing the lattice potentials. If the real 
electron and photon are generated in the way described 
in this paper, the SUSY breaking might be related to the 
change of space-time structure in early universe, which 
has been beyond the present model. 

We consider an integer boson filling n a —1 for Ubb ^ 
I. In this region, since n^ 1 and n" 5 ^ 1 are very small, 
there are no solutions of the mean field equations with 
ya^i -> q f or ^ e cr jti ca l temperature equations. Thus, 
all slave bosons and CFs with a ^ 1 are confined. Only 
T 1 > can be found. Below T 1 , the CF c x and the 
slave boson h\ are deconfined in the mean field sense. 
We shall see the CF C\ may be identified as the 'elec- 
tron' while hi as the Higgs boson. In the inset of Fig. 
1, we plot for a set of given parameters. The curve 
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T p {Ubb) is corresponding to /j,p(T p ) = 0, i.e., the effec- 
tive chemical potential of c\ vanishes. For the fermion 
half filling, f-ip(T = 0) — ► 0. This means that the Fermi 
surface of CF c\ shrinks and implies an insulating state, 
a checkerboard crystal as we shall see soon. Meanwhile, 
we also find — for the half-filling of the fermions. 

We now go to concrete mean field solutions and focus 
on the near half filling of the fermions. There is no flux 
mean field state in the present model because x h ' c are 
real. Two possible solutions are the dimer and uniform 
phase. The dimer phase may be favorite in the fermion 
half-filling. However, sightly away from the half filling, 
the most stable mean field state is the uniform phase, 
as that for the t-J model H, H^ . say, for a = 1, each 
bond is endowed with the same real values of Xa h if a 
CF (i.e., a\fj\0) = c^|vac)) is surrounded by six bosons 
(i. e., a]|0) = /ijjvac)). In the fermion half filling, this is 
a checkerboard crystal state 0. Slightly away from the 
half filling, this is a UMF state with h\-c\ bond, similar 
to that in the fermion Hubbard model or t-J model (see, 
e.g., 0). In the uniform phase, we choose Xa > 0- The 
dispersions of the CFs and slave bosons are £a' h (k) = 
T^fIx^' ! Si cos A:,. The CF ground state is around k = 
while the slave boson's is around k = if. 

To see the stability of the UMF state, one has to add 
back the gauge fluctuations on the mean field solution. 
The zeroth component Am restores the exact constraint 
of one slave particle per site. This fluctuation may un- 
stabilize the mean field state. Before integrating over h 
or c, the system is in the strong coupling limit because 
the gauge field has no dynamics, i.e., the Lagrangian of 
the gauge field is L g = — lim^oo -F^F^ (in the con- 
tinuous limit). This means even the slave particles have 
a finite mass in T < T C Q , the system is still in charge con- 
finement phase. Near k — 0, one has to integrate away 
h\. Then, the coupling constant of the gauge field be- 
comes finite. In the continuous limit, the effective gauge 
action reads 



S[A] 



^ 2 

n 



- / d 3 kA^k,u n )A„(k,u n )TlZ(k,Lj n ), (7) 



where is the response function integrating over hi . 

The long wave length behaviors of these response func- 
tions are strongly dependent on the phase structure of the 
constitution bosons. Introduce a Hubbard-Stratonovich 
field $i = V y/a + 1 [c^-i-n ., : c a ,i + /i* ^ A h a ,i] to decouple 
t b terms jjj, E3 • ma y be thought as the order param- 
eter field of the Bose condensation. The phase boundary 
of the boson superfluid (BSF)/normal liquid phases can 
be determined by G^^O) =00 where the $-field 
propagator is given by 



G 1 (k, iu n ) = £ k -el^2{a + l] 



nr — n 



a + l 



fiB — olUbb — Ubf 



iu n + [IB - CkUbB 

(8) 
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FIG. 1: The phase diagram of the CFs in the Ubb-T plane 
for n a = 1, Ubf = Ubb/2 i F = n f = 0.55. The solid 

curves are the critical temperature T** after considering the 
gauge fluctuation and the shrink of the CF Fermi surface. 
The dash lines and the slashed area is the crossover from the 
'electron' gas to the CF domain. The inset is the mean field 
phase diagram for both bosons and fermions. 



As expected, the phase diagram of the constitution boson 
consists of the Bose superfluid (BSF), the normal liquid 
and the MI, in which the MI phase only exists in the 
zero temperature and an integer boson filling factor . In 
the inset of Fig. 1, we show an example for the same 
parameters as the CF mean field phase diagram. 

We now go back to the response function n^„. The 
current-current response function has the form Ilfj — 
(Sij — kikj/k 2 )HL + (kikj /k 2 )"n.T- In the normal liquid 
and MI of the boson, the longitudinal part Hl{0, 0) = 
and the gauge field is massless. If hi condenses, 
U L (k,w) -> K/i 1 )! 2 + when fe,u -> 0. This gives the 
gauge field a mass. This is well-known Higgs mechanism. 
In this sense, hi may be identified as the Higgs boson. 
This happens in the BSF. 

The density-density response function n^} '' of free 
slave boson yields an effective attraction between CFs 
in the long wave length limit. This means that the UMF 
state is not stable if ts term does not contribute to n^ , 
which is indeed the case in the BSF phase. In the MI, the 
ts term may be thought as a perturbation if Ubb and 
\Ubb — Ubf \ are much larger than tg. For the nf = 1 MI 
phase, the second order perturbation and the MI char- 
acter, + = 1, ts terms contribute an effective 
interaction between the CFs or slave bosons 



J 



Jijnfn^ + const 



Jijn^n) + const', (9) 



with Jij = jj—qj 



T 

BB 



We see that this interaction 



is attractive if Ub f > Ub b while it is repulsive if Ubf < 



Ubb- For the repulsive one, if ts , 



for a special spatial structure, the effective interaction 
between CF or the slave boson in the continuous limit is 



4 



a pure Coulomb one J = J2k,k',q ^(9) c k+q c k c k'- q c k' = 
Y,^k', q V (l) h l+c l h ^ h l-c l h ^ with 

47TP 2 

V(q) = — , (10) 

and the 'charge' of the CF is defined by e 2 = 
77 — Try? — ^775 — 7. This repulsive interaction keeps the sta- 

Ubb(U bb -U bf ) r f 

bility of the UMF state against the gauge fluctuations. 
Eq. I|1U|) means that the CFs in the MI phase behave like 
'electrons'. 'Electrons' are a Coulomb gas. Notice that 
the electron gas at the half filling (n? = 1/2) turns to a 
checkerboard 'electron' lattice. In this way, we see that 
the 'electron' and 'photon' (the gauge field) emerge from 
a neutral SUSY particle model. 

In the ulrtacold Bose-Fermi mixture, instead of a pure 
Coulomb interaction caused by a long range hopping 
ts,ij, the interaction between nearest neighbor sites dom- 
inates which is caused by ts,ij = ts =/= if only is 
a nearest neighbor pair. The interaction between CFs at 

T = is given by V(q) = u^0%p y (Also see 0.) 

According to these discussions, we now may figure out 
the phase diagram of the CF in Fig. 1. The mean field 
phase transition temperature is suppressed greatly to 
Tj}* which is determined by the T c and T p . The dash line 
is the estimated crossover line from the 'electron' gas (the 
MI of bosons) to CF domains which arises from the ef- 
fective attraction between CFs. In the BSF, the induced 

attraction between CFs, (n^ (0) (q, 0))" 1 r -^ff-, may 

lead to a p wave superconducting ground state. However, 
this might be difficult to reach in the ultracold mixture 
0. 

The experimental implications of the UMF phase are 
discussed as follows. We consider the 'electron' response 
to an external 'electric' field, 'made' by a change of the 
lattice potential of the fermions. This disturbs the den- 
sity of fermions with H'(t) = — JY wfa, t)(p(ri, t). In a 
time flying experiment, the difference between disturbed 
and undisturbed fermion densities by external field is 
given by n/(r) — n°(r) = (^-)\wf(k = ^y^)| 2 i5n/(k = 



) where t is the flying time, Wf is the Fourier com- 
ponent of the fermion Wannier function and (Sny(k) oc 
n o(k, 0) . If T > Tl* , the density response of the system 
is simply given by the free fermion one, n£, (0) . For T < 
since n* f« n*, the two CF response function is given 

bylE^ 1 = (n^^MlCV 1 - Negative < (0) leads to 
the instability of the CF. However, near the MI, the RPA 
gives n^ 1 = (nS 0) ) - 1 + ttl^ ' ) - 1 + V(q) . The repulsive 
between CFs stabilizes the CF against the gauge fluctua- 
tion. A better experimentally measurable quantity is the 

visibility V = ^(r™)+n/(rll°j wnere r max and r min are 
chosen such that the Wannier envelop is cancelled [Til] . 
The difference between the disturbed and undisturbed 
visibility may directly correspond to the response func- 
tion because n/(r max ) + n/(r min ) w n° f (r max ) + n a f (r min ) 
in denominator. 

The Fermi surface of pure cold fermions has been ob- 
served in a recent experiment Q by the time flying ex- 
periment. We suggest to do the same observation to the 
mixture. It is expected that instead of the constitution 
fermion's Fermi surface, one may observe the 'electron' 
Fermi surface in the 'electron' gas. 

In conclusions, we deduced the low energy physics of 
the SUSY Hubbard model. The 'electron' and 'photon', 
as well as the Higgs boson, were thought as emergent 
objects. The phase diagram was depicted, which showed 
what happened as the temperature was from high to low. 
Possible experiments to verify this theory were suggested 
by means of the cold Bose-Fermi atom mixture. We ex- 
pect to generalize our theory to include a non-abelian 
gauge field |l3l | coupled to fermions with spin degrees of 
freedom. However, it will not be a simple generaliza- 
tion of the U(l) theory. Many problems, such as quark 
asymptotic free in an ultraviolet limit and confinement 
at low energy, have to be solved. The origin of the chi- 
ral fermions in SU(2) x U(l) weak-electric interaction is 
also non-trivial. We expect a relativistic version of the 
present theory. 
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